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ABSTRACT 



A FORTRAN COntrol Program for Engineering Synthesis 
(COPES) has been developed for solving engineering design 
problems. The program maximizes or minimizes a numerically 
defined objective function subject to a set of inequality 
constraint functions. COPES uses the optimization program, 
CONMIN, which includes the conjugate direction method of 
Fletcher and Reeves for unconstrained function minimization 
and a modification of Zoutendijk's method of feasible direc- 
tions for constrained function minimization. Additionally, 
approximation techniques are available for use in optimiza- 
tion, and trade-off studies may be performed. A simple 
design example demonstrates the program capabilities. 
Programming guidelines are presented followed by sample 
input data and output for each program option. 
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I. INTRODUCTION 



Most design processes require the minimization or maximi- 
zation of some parameter which may be called the design ob- 
jective. For the design to be acceptable, it must satisfy a 
variety of physical, aesthetic, economic and, on occasion, 
political limitations which are referred to here as design 
constraints. While part of the design problem may not be 
easily quantified, most of the design criteria can be des- 
cribed in numerical terms. 

To the extend that the problem can be stated in numerical 
terms, a computer program can be written to perform the 
necessary calculations. For this reason, computer analysis 
is commonplace in most engineering organizations. For example, 
in structural design, the configuration, materials and loads 
may be defined and a finite element analysis computer code 
is used to calculate stresses, deflections and other response 
quantities of interest. If any of these parameters are not 
within the prescribed bounds, the engineer may change the 
structural member sizes and re-run the program. The computer 
code therefore provides only the analysis of a proposed design, 
with the engineer making the actual design decisions. This 
approach to design, which may be called computer-aided design, 
is commonly used today. 



6 



Another common use of analysis codes is in trade-off 
studies. For example, an aircraft trajectory analysis code 
may be run repetitively for several payloads, calculating the 
aircraft range, to determine the range-payload sensitivity. 

A logical extension to computer-aided design is fully 
automated design, where the computer makes the actual design 
decisions, or to perform trade-off studies with a minimum of 
man-machine interaction. The purpose of the COPES program 
is to provide this automated design and trade-off capability. 
The user must provide a FORTRAN analysis program for analysis 
of the particular problem being considered. This analysis 
program is written according to a simple set of guidelines 
so that it can be easily coupled to the COPES program for 
automated design synthesis. 

This document describes the capabilities of the COPES 
program and its usage. A simple design example is first pre- 
sented to demonstrate the program capabilities. Guidelines 
are given for writing analysis codes which can be coupled 
directly to COPES. Finally the data organization is outlined 
and sample data is presented. 

The most commonly used option of the COPES program will 
be for design optimization. Two approaches are available 
for this purpose. The first (NCALC = 2) is direct optimiza- 
tion of the function by the CONMIN optimization sub-program 
[Ref. 1] . An alternative to this is through the use of 
approximation techniques (NCALC = 6, Ref. 2). This second 
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option is usually more efficient for problems of under 6 
design variables but which requires costly analysis, espe- 
cially when multiple optimizations are to be performed. 



DESIGN EXAMPLE 



Assume it is required to 
design the cantilevered beam 
shown in Fig. 1. The objec- 
tive is to find the minimum 
volume of material which will 
support the concentrated load. 
That is; 



P=10,000 lb. 




| ** — L=200 i n. — ► | 
E=30x1 0 6 psi 



□ 



H 



Minimize Volume = B‘H*L 



( 1 ) 



Figure 1 



Cantilevered 

Beam 



The bending stress in the beam must not exceed 20,000 psi; 



a 



b 



Me 

I 



6PL 



2 < 



BH - 



20,000 



( 2 ) 



The shear stress must not exceed 10,000 psi; 



v 



3 P 
2 A 




10,000 



(3) 



and the deflection under the load must not exceed one inch; 



6 



PL 3 

3EI 



4PL' 



EBH 3 - 



< 1.0 



(4) 



Additionally, geomet 
so that; 

0.5 £ B 
1.0 < H 
H/B 



ic limits are impo 

£ 5.0 
£ 20.0 
< 10.0 



ed on the beam size 

(5) 

( 6 ) 

(7) 
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Equation (1), the design objective, provides a measure of 
the efficiency of the design, while Eqs . (2-7) define the 
criteria which the structure must satisfy, and are referred 
to as constraints. 

This design problem may be stated in general form; 

Minimize F(X) 

Subject to: 

Gj(X) < 0 i = 1 ,m 

X* < X. < x u 
1 — 1 — 1 

where X is a vector containing the design variables B and H 
and Gj (X) are the constraints defined by Equations (2-7). 
There are eight constraints on the design. The objective 
and constraints are functions of the design variables. 

While this problem is straightforward, its solution is 
not trivial because it is not known which constraints will 
be critical (i.e. which Gj(X) = 0) at the optimum. 

One design approach is to assume that the bending and 
the displacement constraints are critical, so that Equations 
(2) and (4) are equalities, and then solve the two simulta- 
neous equations for the design variables B and H. The 
remaining constraints are then checked to be sure they are 
satisfied, and if not, they are used to obtain a new design 
satisfying all constraints. When a design is obtained where 
two Gj (X) = 0 and the remaining Gj (X) £0, it is an accept- 
able design. 
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The design obtained here yields a volume of 6750 and 
satisfies all constraints. Note, however, that the objec- 
tive function of Equation (1) played no part in the choice 
of design variables so that, using this approach, there is 
no assurance that a minimum volume design will be produced 
[Ref. 3] . More importantly, it is desirable to devise 
techniques for optimum design of systems which may be de- 
fined by more than two variables and by much more complex 
analysis, where the optimum design cannot be determined by 
inspection. 

The COPES program provides this general capability by 
the use of the optimization program CONMIN [Ref. 1] . To use 
this design capability, a FORTRAN code must be provided 
which will calculate the various parameters. In writing the 
analysis code; 1) it is written in subroutine form with 
SUBROUTINE ANALIZ (ICALC) as the main routine, 2) it is seg- 
mented into INPUT, EXECUTION AND OUTPUT and 3) all parameters 
which may be design variables, object functions or con- 
straints are contained in a single labeled common block called 
GLOBCM. 

To demonstrate the simplicity with which a design- 
oriented analysis code can be written, the following FORTRAN 
subroutine was produced for the analysis of the cantilevered 
beam in Figure 1. 
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SUBROUTINE ANAL I Z (ICALC) 

COMMON /GLOBCM/ B ,H , VOL , BSTRES , SHRSTR, DELTA ,HB ,E ,AL 
IF (ICALC.GT.l) GO TO 10 
C --- INPUT OR INITIALIZATION. 

B= 2 . 5 
H=10 . 

P=10,000. 

E=30 .E+6 
AL=200 . 

WRITE (6,30)AL,P,E ,B,H 
RETURN 

C --- EXECUTION. 

10 CONTINUE 

VOL=AL*B*H 

BSTRES=6 .*P*AL/(B*H*H) 

SHRSTR=1 . 5*P/ (B*H) 

DELTA=4 .*P*(AL**3) / (E*B*(H**3)) 

HB=H/B 

IF (ICALC.LT. 3) RETURN 
C --- PRINT RESULTS. 

20 WRITE (6 , 30) AL , P ,E ,B ,H 

WRITE (6,40) VOL , BSTRES , SHRSTR, DELTA ,HB 
30 FORMAT(/////5X, 17HCANT I LEVERED BEAM//5X, 8HAL =,F9.3/5X, 

* 8HP =,E12.5/5X, 8HE =,E12.5//5X, 8HB =,F9.3/5X, 

* 8HH =,F9.3) 

40 FO RMAT ( / 5 X , 8HVOL = , F9 . 3//5X , 8HBSTRES= , E 1 2 . 5/5X , 

8HSHRSTR= 

* E 12 . 5/5X , 8HDELTA= ,E12 . 5/5X , 8HH/ B =,F9.3) 

RETURN 

END 



This routine may be executed as a simple analysis program 
by the following main program; 

C MAIN PROGRAM TO EXECUTE SUBROUTINE ANALIZ. 

DO 10 I = 1,3 
10 CALL ANALIZ (I) 

STOP 

END 



Moreover, ANALIZ can be coupled directly to the COPES 



program to perform this same function, or to perform optimiza- 
tion or trade-off studies. 



This subroutine was coupled to COPES and an optimization 
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was performed (NCALC = 2) to yield the following design; 



CANTILEVERED BEAM 



L 

P 

E 

B 

H 



200.00 



. 10000E + 05 
. 30000E+08 
1.82 



VOL 

BSTRES 

SHRSTR 

DELTA 

H/B 



18.16 

6607.61 



. 20002E + 05 
. 45402E + 03 
. 97906E+00 
9.98 



The parameters L, P, E, B and H were input to the pro- 
gram and the remaining parameters were calculated. The design 
variables B and H were changed during the optimization process 
to obtain this design. 

This design was achieved with 45 calls to ANALIZ with 
ICALC = 2 (45 analyses). This design could surely have been 
found with fewer analyses were it performed by hand calcula- 
tions. However, once having written the analysis subroutine, 
numerous other designs may be obtained for different materials, 
loading, or design stresses with only minimal effort. 
Furthermore, the design obtained here of volume = 6607.6 is 
very near the theoretical optimum of 6603.9, while that which 
satisfied stress and displacement constraints simultaneously 
was not. 

This very simple design example underscores the power of 
numerical optimization techniques and the ease with which they 
may be applied. The key to efficient use of the COPES program 
is the requirement that the ANALIZ code be written in a standard 
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format so that it may be coupled to COPES without modifica- 
tion. The following section contains guidelines for writing 
a design-oriented analysis code. 
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II. PROGRAMMING GUIDELINES 



In developing any computer code for engineering analysis, 
it is prudent to write the code in such a way that it is 
easily coupled to a general synthesis program such as COPES. 
Therefore, a general programming practice is outlined here 
which in no way inhibits the use of the computer program in 
its traditional role as an analytic tool, but allows for sim- 
ple adaptation to COPES. This approach is considered good 
programming practice and provides considerable flexibility 
of design options. Only five basic rules must be followed: 

I. Write the code in subroutine form with the primary 
routine called as; SUBROUTINE ANALIZ (ICALC) The 
name ANALIZ is compatable with the COPES program and 
ICALC is a calculation control. Note that subroutine 
ANALIZ may call numerous other subroutines as re- 
quired to perform the necessary calculations. 

II. Segment the program into INPUT, EXECUTION and OUTPUT. 
The calculation control, ICALC, will determine the 
portion of the analysis code to be executed. 

ICALC = 1; the program reads all data required to perform 
the analysis. Also, any initialization of constants which 
will be used repetitively during execution is done here. 

This initial input information is printed here for later re- 
ference and for program debugging. 
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ICALC = 2: the program performs the execution phase of 

the analysis task. No data reading or printing is done here, 
except on user-defined scratch disc. Data may be printed 
here during program debugging, in which case it should be 
controlled by a print control parameter which is read during 
input. In this way, this print may be turned off after the 
program is debugged, but may be used again during future 
program expansion and debugging. The reason that printing 
is not normally allowed during execution is that when 
optimization is being done, the code will be called many 
times with ICALC = 2, resulting in voluminous print. 

ICALC = 3: the results of the analysis are printed. 

Also the essential input parameters which may have been 
changed during optimization should be printed here for easy 
reference. In some cases, so much information is generated 
when ICALC = 2 that it would be inefficient to store it 
internally or on disc for printing when ICALC =3. In this 
case, it may be desirable to actually execute when ICALC = 3 
with a print code turned on to print his intermediate infor- 
mation. However, this approach should be avoided because 
it requires an additional complete execution of the program. 

In summary, when; 

ICALC = 1 Read input data. 

ICALC = 2 Execute the analysis. 

ICALC = 3 Print the results. 
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III. Store all parameters which may be design variables, 
objective functions or constraints in a single 
labeled common block called GLOBCM. The order in 
which they are stored is arbitrary. A listing of 
the COPES program should be checked to see how many 
parameters my be stored in GLOBCM (the dimension of 
ARRAY) . Initial distribution of COPES allows for 
1500 parameters. 

IV. During execution or output, no parameters which are 
read during input should be updated. For example, 
if variable X is initialized during input, the 
execution segment must not update X such as X = X + 
3.2. Instead a new variable, Y = X + 3.2 should be 
defined. 

V. Write all programs in standard language, avoiding 
machine dependent capabilities such as multiple 
entry point (IBM) , DEFINE statements (UNIVAC) and 
seven letter FORTRAN names (CDC) . While this guide- 
line is not essential to the use of the analysis 
code within the COPES program, it makes the analysis 
code much more transportable between different 
computer systems, a capability which easily justi- 
fies a slight reduction in efficiency on a given 
machine . 

Adherence to these guidelines not only leads to a more 
readable and machine independent computer code, but allows 
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this code to be coupled to the COPES program without 
modification . 

Having written the analysis code, it may be executed 
either with a simple main program or within the COPES program 
to perform the analysis. To insure that guideline IV is 
followed, the following main test program is recommended. 

Note that this program calls ANNALIZ twice with ICALC = 2 and 



ICALC = 3, to show that the same result is obtained repetitively 

C MAIN PROGRAM TO CHECK SUBROUTINE ANALIZ , 

C READ, EXECUTE AND PRINT. 

DO 10 ICALC =1,3 
10 CALL A.NAL I Z (ICALC) 

C EXECUTE AND PRINT AGAIN TO BE SURE THE RESULTS 
C DO NOT CHANGE. 

DO 20 ICALC = 2, 3 
20 CALL ANALIZ (ICALC) 

STOP 

END 

This program was executed with the ANALIZ subroutine for 
the cantilevered beam example to yield the following result: 



CANTILEVERED BEAM 



AL 



P 



H 



E 



B 



200.00 
. 10000E+05 
. 30000E + 08 
2.00 
5.00 
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CANTILEVERED BEAM 



AL = 200.00 

P = . 10000E+05 

E = . 30000E+08 

B =2.00 

H =5.00 

VOL = 2000.00 

BSTRES = . 24000E+06 

SHRSTR = . 15000E+04 

DELTA = . 42667E+02 

H/B =2.50 

CANTILEVERED BEAM 



AL 

P 

E 

B 

H 

VOL 

BSTRES 

SHRSTR 

DELTA 

H/B 



= 200.00 
= . 10000E+05 
= . 30000E+08 
= 2.00 
= 5.00 
= 2000.00 
= . 24000E + 06 
= . 15000E+04 
= . 42667E+02 
= 2.50 



This design was used as the initial design for the optimi- 
zation presented previously. Note that, while the volume here 
is lower than the optimum, the bending stress and displacement 
each exceed the imposed limits by more than an order of 
magnitude . 

Once the analysis code has been written, it can be coupled 
to the COPES program without modification. If it is desired 
to perform a simple analysis using COPES, only three data 
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cards are required for the COPES program, namely a TITLE card, 
a control parameter, NCALC = 1, and an END card. If the 
optimization or parametric analysis (sensitivity) capabilities 
of COPES are to be used, additional data must be read. This 
data will identify which parameters in the global common block, 
GLOBCM, are used. To set up the COPES data, the user must 
have a basic understanding of how the data in the global com- 
mon block is accessed by COPES. This is outlined in the 
following section. 
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III. DATA MANAGEMENT 



In order to perform design operations, the COPES program 
must access the data in common block GLOBCM. This is done by 
defining the location in GLOBCM where a specified parameter 
resides. For example, consider the common block for the can- 
tilevered beam; 

COMMON/ GLOB CM/ B ,H, VOL ,BSTRES ,SHRSTR, DELTA, HB ,E ,AL 

The volume of material, VOL, is the third parameter in 
the common block; that is, it resides in location 3, referred 
to as the global location number. Similarly the bending 
stress, BSTRES, is in global location 4 and the beam width is 
in global location 1. Thus, the parameters are referred to 
by their respective location numbers in global common. 

For convenience in preparing data for the COPES program, 
a simple "CATALOG" of parameters may be defined. For the 
cantilevered beam, this catalog would be; 

GLOBAL FORTRAN 

LOCATION NAME DEFINITION 



1 


B 


Beam width 


2 


H 


Beam height 


3 


VOL 


Volume of Material 


4 


BSTRES 


Maximum bending stress 


5 


SHRSTR 


Maximum shear stress 


6 


DELTA 


Deflection under the load 


7 


HB 


Ratio, H/B 
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GLOBAL 


FORTRAN 




LOCATION 


NAME 


DEFINITION 


8 


E 


Young's modulus 


9 


AL 


Length of beam 



As another example, consider a global common block con- 
taining arrays; 

COMMON/ GLOBCM/A, Y(10), Q, C(2,2), H 

The variable catalog for this common block is; 



GLOBAL 

LOCATION 


FORTRAN 

NAME 


DEFINITION 


1 


A 


Area 


2 


Y ( 10) 


Vector of y-c 


12 


Q 


• 


13 


C(2 ,2) 


• 


17 


H 


etc . 



The dimensions are given with the FORTRAN name as a 
reminder that the parameter is an array. In this case, the 
third parameter in the Y array is in global location 4. 
Remembering that arrays are stored column by column the 
C(l,2) array location is in global location 15. 

It will be seen that identifying parameters according 
to their location in GLOBCM provides a great deal of flexi- 
bility in using the COPES program for design. 

Appendix A contains blank forms for writing the variable 
catalog for the user's particular problem. 
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IV. COPES TERMINOLOGY 



The copes program currently provides six specific 
capabilities ; 

1. Simple analysis, just as if COPES was not used. 

2. Optimization - Minimization or maximization of one 
calculated function with limits imposed on other 
functions using the CONMIN subprogram. 

3. Sensitivity analysis - the effect of changing one 
or more design variables on one or more calculated 
functions . 

4. Two- variable function space - analysis for all 
specified combinations of two design variables. 

5. Optimum sensitivity - same as sensitivity analysis 
except at each step, the design optimized with 
respect to the remaining independent design vari- 
ables. 

6. Approximate optimization - optimization using 
approximation techniques. Usually more efficient 
than standard optimization for up to 6 design 
variables or if multiple optimizations are to be 
performed . 

In defining the data required to execute the COPES pro- 
gram, the following definitions are useful: 

Design Variables - Those parameters which the optimiza- 
tion program is allowed to change in order to improve the 
design. Design variables appear only on the right hand side 
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of equations in the analysis program. COPES considers two 
types of design variables, independent and dependent. If 
two or more variables are always required to have the same 
value or be in a constant ratio, one is the independent 
variable while the remaining are dependent variables. For 
example, if the height is required to be 10 times the width 
of the cantilever beam, B would be the independent variable 
while H would be the dependent variable. 

Objective Function - The parameter which is to be mini- 
mized or maximized during optimization. Also the parameters 
calculated as functions of specified design variables during 
a sensitivity or two-variable function space study. Objective 
functions always occur on the left side of equations, unless 
the objective function is also a design variable (the beam 
height may be minimized as an objective function if it is 
also a design variable. In this case the minimum height is 
found for which no constraints are violated). An objective 
function may be linear or non-linear, implicit or explicit, 
but must be a function of the design variables to be meaning- 
ful. 

Constraint - Any parameter which must not exceed speci- 
fied bounds for the design to be acceptable. Constraint 
functions, always appear on the left side of equations. Just 
as for objective functions, constraints may be linear or 
non-linear, implicit or explicit, but must be functions of 
the design variables. 
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Constraint Set - A group of constraints which appear 



consecutively in the global common block and which all have 
the same limits imposed. This is a convenience which allows 
several constraints to be identified with a minimum of data. 

Global Common - Common block GLOBCM containing design 
information. 

Global Location - Location of a particular parameter in 
GLOBCM. 
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V. COPES DATA 



The COPES program reads data from unit S and writes output 
on unit 6. Units 20 and 40 are used as scratch files. The 
scratch file numbers may be changed by changing two cards at 
the beginning of the COPES program. 

In order to execute the COPES program it is necessary to 
provide formatted data for COPES, followed by data for the 
ANALIZ program which is coupled to COPES. This section de- 
fines the data which is required by COPES. The data is seg- 
mented into "BLOCKS" for convenience. All formats are alpha- 
numeric for TITLE, and END cards, F10 for real data and 110 
for integer data. 

The COPES data begins with a TITLE card and ends with an 
END card. This is followed by data to be read by the user 
supplied subroutine ANALIZ or when ICALC = 1. 

Comment cards may be inserted anywhere in the COPES data 
stack prior to the END card, and are identified by a dollar 
sign ($) in column 1. 

Data coding forms are provided in Appendix B. 
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VI. UNFORMATTED DATA INPUT 



While the user’s sheet defines COPES data in formatted 
fields of ten, the data may actually be read in a simplified 
fashion by separating data by commas or one or more blanks. 

If more than one number is contained on an unformatted data 
card, a comma must appear somewhere on the card. If expo- 
nential numbers such as 2.5+10 are read on an unformatted card, 
there must be no embedded blanks. Unformatted cards may be 
intermingled with formatted cards. Real numbers on an un- 
formatted card should have a decimal point. 

Examples : 

Unformatted data; 

5,7,1. 3,1. 0+20,0, -5.1 
5,7,1.3,1.0+20, ,-5.1 

571.3 1.0+20, , -5.1 

5 7 1.3, 1.0+20 0 -5.1 

Equivalent formatted data; 

col+ 10 20 30 40 50 60 70 80 

5 7 1.3 1.0+20 0 -5.1 

Unformatted data; 

2 

2.3 
2 3 
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Equivalent formatted data; 

col-*- 10 20 30 

2 

Note: This data has been right 

2 3 justified. 

2 3 Note: This data contains no commas, 

so it is assumed to be format- 
ted already. 

Unformatted data; 

1.2.3.4.5.6.7.8.9.10.11 
Equivalent formatted data; 

col-*- 10 20 30 40 50 60 70 8_0 

12345678 

9 10 11 

Note that two formatted data cards are created here. 

Unformatted data; 

1,2 , 3, 4, 5, 6 , 

7.8.9.10.11 

Equivalent formatted data; 

col+ 10 20 30 40 50 60 60 80 

1 2 3 4 5 6 

7 8 9 10 11 

Note that the above two examples do not produce the same for- 
matted data cards. 
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DESCRIPTION: Title card. 
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Two variable function space. Data of blocks A, B, and R-V is required. 
Remaining data is optional. 
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Print above plus gradient information. 

Print above plus each proposed design vector, objective function 
and constraint values during the one -dimensional seaich. 



ITMAX: Maximum number of optimization iterations allowed. DEFAULT 
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DESCRIPTION : Floating point optimization program parameters. 
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FIELD CONTENTS 

3 CT : Constraint thickness parameter. DEFAULT = -0.1. 

4 CTMIN: Minimum absolute value of CT considered in the optimization 

process. DEFAULT = 0.004. 
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1) The DEFAULT values for these parameters usually work well. 
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DATA BLOCK H OMIT IF NDV = 0 IN BLOCK 
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If two or more adjacent parameters in the global common block have the same 
limits imposed, these are part of the same constraint set. 
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DATA BLOCK K OMIT IF NXAPRX = 0 IN BLOCK 
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XFACT1: Multiplier on DXI when the diagonal elements of the H matrix are 

available. Default = 1.5. 

XFACT2 : Multiplier on DXI when all elements of the H matrix are available. 



DATA BLOCK L OMIT IF NXAPRX = 0 IN BLOCK B OR IN X LOC = 0 IN BLOCK _J 





H 


O 






3 


rH 
k— 1 






2 


OO 






o 








Ph 










. 








• 






00 


• 








• 








• 








• 

• 

• 




CO 








CD 








rH 
















oj 








*H 








U 


vO 






a3 








> 








SO 








C 








•H 








+-> 








aj 








6 


LO 






•H 








X 








o 








5h 








Pu 








Ph 








a3 




''H- 








X 




HH 




CJ 




O 




o 








kJ 




to 








C 








o 








•H 








4-> 




to 




o3 




X 




O 


to 


CJ 




O 




O 




rH 








rH 








03 


• • 






03 


P 


PJ 








O 






03 




rH 


Pu 




X 




cj 


l-il 


03 


CJ 




< 




o 




X 








• • 


W 


















o 


Q 






rH 


k— 1 










H 


< 




rH 




PU 






X 




k— 1 


H 


rH 


CJ 




X 






o 




CJ 






t-J 




CO 


pci 








LU 


O 








Q 1 Hh 








CD 

rH 

XI 

a 

•H 

5h 

o3 

> 

SO 

c 

•H 

+-> 

o3 

£ 

•H 

X 

O 

Sh 

Cu 

Pu 

o3 







<4H 

O 



c 

o 



oj 

o 

o 



o3 

O 



cj 



X 

CJ 

o 

hJ 



°0 



I 




CO 

X 

X 




44 



If INXLOC = 0 in BLOCK J, this data is not read. In this case, the data is 
defaulted to be the global locations of the design variables ( IDSGN values 
in BLOCK G) . 



DATA BLOCK M OMIT IF NXAPRX = 0 IN BLOCK B OR INFLOC = 0 IN BLOCK 
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If INFLOC = 0 in BLOCK J, this data is not read. In this case, the data is 
defaulted to be the global locations of the objective function (IOBJ in BLOCK 
E) followed by the global locations of the constrained parameters (ICON-JCON 
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DATA BLOCK 0 OMIT IF NXAPRX = 0 IN BLOCK B OR NXFS = 0 IN BLOCK 
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to contain data. 
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REMARKS 

1) More than eight values of the sensitivity variable are allowed. Add data 
cards as required to contain the data . 



DATA BLOCK R OMIT IF N2VAR = 0 IN BLOCK 



PQ 



I 



< n 
u 
o 
+-> 
o 
e 

cd 

cd 

a 



o 

Jh 

4-> 

C 

O 

u 

CD 

u 

cd 

p. 



to 



C 

o 



L) 

rt 

D 

m 

(D 

rH 

X 

cd 

•H 

cd 

> 

O 

£ 

H 



x 

o 

h— I 

H 

Pin 

I— I 

PJ 

CJ 

CO 

pq 

Q 



W 

-J 

pq 

pq 

Q 

< 



o 

Uh 




a 






C 




o 






o 




•H 






•H 




+-> 






4 -> 




CJ 






CJ 




a 






c 




D 






D 




m 






m 




CD 






CD 




rH 






rH 




x> 






JQ 




cd 






cd 




•rH 






•H 




Jh 






Jh 




cd 






cd 




> 






> 




o 






o 




£ 






£ 










4 -> 




CD 






CD 




jq 




• 


jq 




+-> 




-d 


+-> 








CD 










Jh 


c 




•rH 




o 


»H 














<D 




•rH 


CD 




rH 




CO 


rH 




X 




jq 


X 




cd 




o 


cd 




♦ H 




o 


•H 




U 






U 




cd 




CD 


cd 




> 




X) 


> 




i 

X 




o 










4 -> 






CD 






CD 




jq 




X 


jq 




4 -> 






4 -> 








m 






m 




o 


Mh 




o 






o 








CO 






c 




CD 


iq 




o 




D 


o 








rH 


•rH 




+-> 




cd 


4 -> 




cd 




> 


cd 




u 






cj 




o 




m 


o 




rH 




o 


rH 




rH 


• 


Jh 


rH 




cd 


CD 


CD 


cd 


CD 


X 


u 


X) 


X 


CJ 


o 


cd 


G 


O 


cd 


rH 


a, 


D 


rH 


P. 


u 


co 




o 


CO 



• • 


• • 


• • 


X 


X 


>* 


> 


> 


> 


CNl 


OJ 


CNI 




S 




rH 


CNI 


K) 



50 



M2VY : Number o£ values of Y to be considered. 

IP2VAR: Print control. If IP2VAR.GT.0, detailed print will be called 

at each step (each X-Y combination). DEFAULT = No print. 
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More than eight functions are allowed. Add data cards as required 
to contain the data. 



DATA BLOCK T OMIT IF N2VAR = 0 IN BLOCK 
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1) More than eight values are allowed. Add data cards as required to contain 
the data. 



X 

nd 

D 

P 

to 

0 

a 

cd 

PD 

(0 

£ 

o 



a 

a 

D 

<4-1 



CD 



X* 

cd 

•H 

P 

cd 

> 

O 

£ 

4-> 



PQ 



cd 



U 

O 

PJ 

P0 

a: 



o 

II 

pc 

<NI 

pp 



H 

i— « 
cei 

o 



c 



CD 

rH 

-Q 

cd 

•H 

P 

cd 

> 

< 

>- 



CD 

SC 



<4-4 

O 

CO 

CD 

3 

r— I 

cd 

> 



13 I 



CJ 

O 

-J 

PQ 



< 

H 

< 

O 



o 

I— I 

H 
Ph 
t — i 

Pi 

u 

CO 

w 

Q 



W 

J 

Ph 

S 

2 

PJ 



H 

I 



O 

PP 



H 


O 




PC 


rH 




O 


PP 




PP 


00 






♦ 




00 


• 






• 










r-. 


• 
















X> 


• 














o 




LO 


• 


LO 


>< 


o 






ro 




H* 


o 


H" 


>* 


• 






vP 






rH 






O 




co 


• 


fO 




CNJ 






rH 




<NI 


o 


<NI 


>-• 


• 






00 






o 


r— 1 


rH 


• 






H* 



£ 



0 

o 

cd 

Pu 

10 

a 

o 

•H 

P 

o 

c 

D 

<4-1 



0 



rQ 

cd 

•H 

P 

cd 

> 

o 

£ 




cd • 
•H 0 
P P 
cd 0 
> Jd 



r d 

cd 

0 0 
rC P 
P 

0 

<4-1 P 

o cd 



PS 

• H 
cd 
p 
PJ 
o 
o 

o 

p 

nd 

0 

P 

• r4 

D 

cr 

0 

p 

CO 

cd 

CO 

nd 

P 

cd 

o 

cd 

p 

cd 

nd 

nd 

nd 

< 



nd 

0 

£ 

o 



cd 

0 

P 

cd 



CO 

0 

D 



cd 

> 



CO CO 


Pd 


0 0 


bO 


D 0 


*H 


rH rH 

cd cd 


0 


> > 


P * 

cd cd 
X *-> 


, . 


P cd 


1— 1 


nd 



>H 0 

P 0 





L0 




« 


Q 


PC 


PI 


CO ^ 


PJ 


i Pi 


i — i 


fH W 


PP 


pc 



53 



DATA BLOCK 
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1) This card MUST appear at the end of the COPES data. 

2) This ends the COPES input data. 

3) Data for the user-supplied routine, ANALIZ, follows this. 



VII. SAMPLE COMPUTER SOLUTIONS 



Sample solutions are presented here for the various COPES 
program options defined by the parameter NCALC in DATA BLOCK 
B. The ANALIZ routine given in the design example section was 
used to produce the results given here. Note that SUBROUTINE 
ANALYZ does not read data, so only COPES data is required here. 
In the usual case where data is read as input to the analysis 
routine, this data would follow directly after the COPES 'END' 
card. 

The output from a COPES program execution includes a title 
page followed by a copy of the input data. Then the required 
program executions are performed and final output information 
is printed. Figures 2-7 contain output for all options of 
COPES as follows: 

PROBLEM SOLVED 



NCALC 

1 

2 

3 

4 

5 

6 



FIGURE 

2 

3 

4 

5 

6 
7 



Analysis only 

Optimization by CONMIN 

Sensitivity study 

Two-variable function space study 

Optimum sensitivity 

Optimization using approximation 
techniques 
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Note that data for the COPES program is read depending on 
the value of the parameters NDV, NSV, N2VAR and NXAPRX in DATA 
BLOCK B. The actual program execution is determined by the 
value of NCALC. Therefore, data may be read for all program 
options even though only one option of the program is to be 
executed. Figure 8 is a copy of the combined data for all the 
previous examples. This data may be used for any program op- 
tion simply by changing the value of NCALC. Figure 9 is a 
copy of this same data using the simplified data input mode 
and Figure 10 is the same data without the comment cards. 
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FIG. 2 - ANALYSIS ONLY 
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CANTILEVERED BEAP 
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FIG. 2 - CONCLD 
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FIG. 3 - OPTIMIZATION 
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IlEVEFED EE AM ANALYSIS AND DESIGN 
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FIG. 3 - CONT 
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FIG. 5 - TWO- VARIABLE FUNCTION SPACE STUDY 
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FIG. 6 - OPTIMUM SENSITIVITY 
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FIG. 7 - OPTIMIZATION USING APPROXIMATION TECHNIQUES 
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ILEVERED EEAK ANALYSIS AND CESIGN 
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FIG. 9 - SIMPLIFIED DATA INPUT 
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FIG. 9 - CONCLUDED 
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SIMPLIFIED DATA INPUT WITHOUT COMMENT CARDS 
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APPENDIX A 
GLOBAL CATALOG FORMS 
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APPENDIX B 
COPES DATA FORMS 
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DATA BLOCK 0 - OMIT IF NXFS 
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